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ABSTRACT 


Parameter and state estimation techniques are discussed for an ellip- 
tic system arising in a developmental model for the antenna surface of the 
Maypole Hoop/ Column antenna. A computational algorithm based on spline 
approximations for the state aud elastic parameters is given and numerical 
results obtained using this algorithm are summarized* 


i. imoeccnon 


Results are presented from a Langley program directed towards 
developing computationally efficient identification techniques for flexible 
systems modeled by partial differential equations with an emphasis on large 
space structures. Initial efforts have been directed towards extending the 
spline-based theory and computational techniques used by the first two 
authors [l]-[6] in solving identification problems with delay and partial 
differential equation models in one spatial variable to solve distributed 
problems in several spatial variables* Additionally, ii. order to support 
Langley’s technology development program [7] in large space antennae, a 
oarameter and state estimation algorithm has been derived for a prototype 
distributed model of the Maypole (Hoop/Column) * antenna reflector 
surface [8]. The next section describes the Hoop/Coluan antenna and pre- 
sents the identification problem being considered. The state and parameter 
estimation approach is then outlined and discussed in tne context of the 
Hoop/Column application. Subsequent sections include mathematical details 
of the antenna application and numerical results. 
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II. THE MAYPOLE (HOOP/OOUMi) AMTEHMA 


For the purpose of technology development, the NASA Large Space 
Systems Technology (LSST) program office has pinpointed focus missions and 
identified future requirements for large space antennas for communications, 
earth sensing, and radio astronomy [7]. In this study, particular emphasis 
is placed on mesh deployable antennas in the 50-120 meter diameter cate- 
gory. One such antenna is the Maypole (Hoop/ Column) antenna shown for the 
100m point-design in Figures 1 and 2. This antenna concept is being devel- 
oped by the Harris Corporation, Melbourne, Florida, under contract to the 
Langley Research Canter [8] . 

The Hoop/ Column antenna consists of a knitted gold-plated molybdenum 
wire reflective mesh stretched over a collapsible hoop that supplies the 
rigidity necessary to maintain a circular outer shape. The annular 
membrane-like reflector surface surrounds a * telescoping mast which provides 
anchoring locations for the mesh center section (Fig. 1). The mast also 
provides anchoring for cables that connect the top end of the mast to the 
outer hoop and the bottom end of the mast to 48 equally spaced radial 
graphite cord truss systems woven through the mesh surface [8] • Tensions 
on the upper (quartz) cables and outer lower (graphite epoxy) cables are 
counter balanced to provide stiffness to the hoop structure. The inner 
lower cables produce, through the truss systems, distributed surface load- 
ing to control the shape of four circular reflective dishes (Figs. 1 and 2) 
on the surface. 

After deployment or after a long period of operation, the reflector 
surface may require adjustment. Optical sensors are to be locat'd on the 
upper mast which measure angles of retroreflective targets placed on the 
truss radial cord edges on the antenna surface. This information can then 
be processed using a ground-based computer to determine a data set of val- 
ues of mesh surface location at selected target points. If necessary, a 
new set of shaping (control) cord tensions can be fed back to the antenna 
for adjustment. 

It is desirable to have an identification procedure which allows one 
to estimate the antenna mesh shape at arbitrary surface points and the 
distributed loading from data set observations- It can also be anticipated 
that environmental stresses and the effects of aging will alter the mesh 
material properties. The Identification procedure must also allow one to 
address this issue* 

Considering the antenna to be fully deployed and in static equilib- 
rium, a distributed mathematical model which describes the antenna surface 
deviation from a curved equilibrium configuration is under investigation 
(for preliminary findings, see [9]). Using a cylindrical coordinate system 
with the z-axis along the mast, it is expected that the resulting model 
will entail a system of coupled second-order linear partial differential 
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equations in two spatial variables. The coefficients of these equations 
are functions of the material properties of the stretched mesh. The deri- 
vation and computer software for this model are still under development. 
In the meantime, a simpler developmental (prototype) problem has been 
solved which is descriptive of the original problem. 

For the developmental problem, the loading is assumed to be normal to 
the horizontal plane containing the hoop rim, and the mesh surface is 
assumed to be described by the static two-dimensional stretched membrane 
equation flO] with variable stiffness (elastic) coefficients and appropri- 
ate boundary conditions for the Hoop/Coluon geometry. Mathematically, in 
polar coordinates, we have 


• 7 It [ rE(r ’ e) It] - 7 le [ E(r > 0) It] • f(r * e) (1) 


where u(r,0) is the vertical displacement of the mesh from the hoop 
plane, f(r,0) is the distributed loading force per unit area, and 
E(r,0) >0 is the distributed stiffness (elastic) coefficient of the mesh 
surface (force/unit length). Equation (1) is to be solved over the annular 
region Q = [e,R] x [0,2*]. Appropriate boundary conditions are 


u(e,0) - u Q 
u(R,0) * 0 

along with the periodicity requirement 

u(r ,0) - u(r ,2*) , 


( 2 ) 


(3) 


where R is the radius from the mast center to the circular outer hoop, 
e is the radius from the mast to the beginning of the mesh surface (see 
Fig. 2), and uq is the coordinate at r » e of the meah surface below 
the outer hoop plane. 

We further assume that the distributed loading along with a data set 

of vertical displacements, u (r ,0 ), at selected points (r ,0 ) on the 

m 1 j J 

mesh surface is known. Given this information, the developmental problem 
is to estimate the material properties of the mesh as represented by 
E(r,0) and produce state estimates of the surface represented by u(r,0) 
at arbitrary (r,0) points within ft. The procedure applied to solve this 
problem is discussed in the next section. 
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III. THE SYSTEM IDENTIFICATION APPROACH 


The first two authors and their colleagues have derived technique? for 
approximating the solutions to systems identification and control probaems 
involving delay equation models and partial differential equation models in 
one spatial variable and have used them in a variety of applications 
{ll] t [12] • The Hoop/Column application requires an extension of the theory 
and numerical algorithms to elliptic distributed systems in several spatial 
variables. The approach, when specialized to the system identification 
problem, may be s t nmarized as follows: (1) select a distributed parameter 
formulation containing unknown parameters for a specific system; (2) mathe- 
matically "project** the formulation down onto a finite dimensional subspace 
through some approximation procedure such as finite differences, finite 
elements, etc.; (3) solve the Identification problem within the finite 
dimensional subspace obtaining an estimate dependent upon the order of the 
approximation embodied in the subspace; (4) successively increase the order 
of the approximation and, in each case, solve the identification problem so 
as to construct a sequence of parameter and state estimates ordered with 
increasing refinement of the approximation scheme; (5) seek a mathematical 
theory which provides conditions under which tue sequence of approximate 
solutions approaches the distributed solution as the subspace dimension 
increases with a convergent underlying sequence of parameter estimates. 

In applying this approach to the developmental problem, the stiffness 
function is parametrized In terms of cubic splines of fixed order; thir 
converts the estimation of E(r,G) into a finite dimensional parameter 
estimation problem. After writing the energy functional generic to the 
membrane equation, the Galerkin procedure is used to project the distri- 
buted formulation onto a finite dimensional state subspace spanned by ten- 
sor products of linear spline functions defined over ft. The approximate 
displacement (state estimate' thus obtained is expressible in terms of the 
spline basis functions. The Galerkin procedure in this case yields alge- 
braic equations which define the displacement approximation coordinates in 
terms of the unknown E(r,8) parameters. In order to solve the approxi- 
mating parameter estimation problem, the parameters defining E(r,0) are 
chosen so that a .east squares measure of the fit error between the 
observed and predicted (by the estimated state) data aet is minimized. 
Finally, following steps (4) aiJ (5) an algorithm is constructed to 
determine the order of the linear spline approximation above which little 
or no further improvement is obtained in the unknown quantities as one 
increases the dimension of the subspaces. Details of this system 
identification approach are presented in the following sections. 
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IV. FINITE DIMENSIONAL APPROXIMATIONS 


Prior to applying the Galerkin procedure [13,14] to perform the finite 
dimensional approximation for the developmental problem, the boundary 
conditions (2) are converted to homogeneous form by introducing the new 
dependent variable 


y(r,e> . u(t,e) - (f-E-§) » 0 


Equation (1) then becomes 

-7 If ('«'•«&)- 7 Is K- e> &) 




e - R 


(4) 


(5) 


with boundary conditions 

y(e,6) - 0 

y(R,6) - 0 (6) 

y(r,0) - y(r,2w) . 


Jollowing the standard formulation (see (13,14]) for the weak or 

A 

variational for- of (5), the energy functional E associated with (5) is 


E(z) - c r & E(r,©)Vz • Vz - (r,0)zj 


rdrd0 


(7) 


where V is the gradient in polar coordinates which, in the form used 
here, li. equivalent to 



1 3_\ T 
r 30/ 


The function f is giver by 


( 8 ) 


?(r,0) 


1 a /rE(r,0)u \ 


(9) 


and the vertical displacement z(r,0) of the mesh surface awa_, from the 
hoop equilibrium plane is a function satisfying the boundary conditions (6) 
and possessing first derivatives on SI in the distributional sense (we 


1 A A 

denote this by ze (ft) = Z). The first variation 6E of E about 

0 , per 

the function y(r,0) is given by 



6E(y;v) 



r,0)Vy 


Vv - ?(r,8)v} rdrd© 


n R 

{E(r,0)Vy • Vr 


[f(r,8)v + E(r ,0)k • Vv] } rdrd0 


( 10 ) 


where 



( 11 ) 


and v is an arbitrary function in Z * H* (0). 

1 O.per 

A 

Given a finite dimensional subspace Z of Z, the Galerkin procedure 

* A 

defines the approximation y as the solution in Z of 



,0)Vy • Vv} rdrdG 


r 2 * 



0)v + E(r,0)k 


Vv} rdrd0 

( 12 ) 


for all veZ. 

For computational efficiency, the basis functions used for the 

a 

representations of y in (12) are taken as tensor products of linear 

A A 

B-splines ([13], p. 27; [14], p. 100). Thus v and y are in the space 
spanned by 

v^ N (r,0) - o”(r)0j(8) , (i - 1,...,M - 1; j - I,...,N) , (13) 

where u^ - a^(r), (i - 1,...,M- 1), and 0^ - 0j(8), ^ ‘ ..,N- 1), 

are standard linear B-spllnes with knots uniformly spaced over fe,R] and 
[0,2*], respectively, modified to satisfy the appropriate boundary 

r Mi 

conditions. The elements [a f are modified to satisfy homogeneous 

N 1 

boun’iry conditions while 0 has been altered to satisfy periodic 

N 

boundary conditions [15]. 
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For y M,N (r,0) 


within the subspace spanned by 



tie can write 


M-l N 


y M,N (r,G) - I l a«(r)w«! N 6j(0) . 
i-1 j-1 J 3 


(14) 


a MM 

Replacing y(r,0) in (12) by y ' (r,0) from (14) and successively set- 
ting v(r,0) ■ v^’ N (r,P) for i * 1,...,M - 1 and j - 1,...,U leads to 

M N 

a set of high-order linear algebraic equations for the w^J coordinates. 


We avoid sparse matrix methods in solving the w 
imposing a separability condition: 


M,N 


U 


equation by 


E(r ,0) - E 1 (r)E 2 (0) . 


(15) 


M M 

As shown in [15], condition (15) reduces the w”’ calculation to the 
so-utio. of the matrix equation ^ 


B M w M ’ N a n + d m w m ’ n c n - e m ’ n 



(16) 

(17) 

( 18 ) 

(19) 

( 20 ) 

(21) 
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and 


-M,N 



f(r,©)a^(r) 0j(9)rdrd9 



c£(r) 


|J rdrdo'j 


( 22 ) 


where, in (17)— (22) , i,p ■ - 1 and j,q » 

Equation (16) is rewritten in the equivalent fora 


[(dVb m ]w m - b ♦ u M ’ N [c H (^)- 1 ] - (d m )- 1 e h 'V 




(23) 


and solved by the Bartels-Stewart algorithm [16]. 

J.n order to estimate, via a numerical scheme, the functional 
coefficients Ej and E 2 , we parametrize these functions so that 

identification is performed over a finite-dimensional parameter set. To 
this end, let 


1 

E. ( r ) * l v l(r) (24) 

k-1 


1 

E (0) - l 6 u (0) (25) 

j-1 * J 


where and 5j are scalar parameters and and y* are 

cubic B-spline functions defined [13, p. 61] over [e,R] and 10, 2*1, 
respectively, whose orders are independent of H and N. The basic spline 
functions are modified so that Mj and its derivatives satisfy periodic 
boundary conditions. 

We turn next to the computer implementation of the identification 
scheme. 


270 



V. COMPUTATIONAL PROCEDURE 


Appealing to the ideas found in previous sections, we now detail an 
algorithm for estimating the coefficients v^, k - 1,..., and 6j, 
j * 1,..., N| , for E(r ,0) that provide the "best fit" between 

estimations of the state u and observed data Uq obtained from various 
sample points on the surface. We may equivalently consider data for y by 
making the transformation 


Vvy ■ wy 



(26) 


for i * 1 ) • • « j L|> find j * 1 9 y Lq • 

A parameter estimation algorithm may je organized into the following 
steps* 


1 • Select an order of approximation for the cubic spline elements 
k * Hi and j ■ 1 Nj , used to 

represent Ej and E 2 * Set n » 1. 


2* Select M and N, number of the linear spline basis elements 
used to represent (and 


3. Assume a nominal set of values for 


v » 



i 


M 


1 


) 


and 


6 - 



(27) 


(28) 


4. Calculate the -oef ficlent matrices in (23) and solve for 
w^»k'(v,6). 

5. Calculate, from (14), y**»N(r.« k ,0j; v,6) and evaluate 

L r L 6 r 

j m, V,<5) - l v,6) - y^vVj * (29) 

6. Proceed to step 8 if jM,N( v> {) is sufficiently small. 

Otherwise, through an optimization procedure, determine a net pair 

* * m N 

(v,6) which decreases the value of J ' . If no such pair can 

be found, go to step 8. 
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7. Set (v,6) ■ (v,6) and return to step 4. 

8. Preserve the current values of and the corresponding 

(v,6) pair as the a*" entry in a sequence of these pairs, 

ordered with increasing H and N. 

9. Proceed to step 10 if sufficient data has been obtained to analyze 
the sequences. Otherwise, set n » n + 1 and return to step 2 
with increased M and N. The current values of (v,6) will be 
used as initial values for the next optimization process. 

10. From analysis of the numerical sequences, select the (M,N) entry 
which indicates the best numerical results. The corresponding 
parameter estimate (v,6) pair yields E(r,6) which determines the 
material properties of the antenna mesh. The matrix 
when used in conjunction with (14), determines a state 
approximation yM,N f or the shape of the antenna surface. 

A convergence theory for the identification algorithm may be found 
in [15]. Numerical results are described in the next section. 


VI. R0W1ICAL RESULTS 


Experimental data for the Hoop/Column antenna is not available at this 
time. Therefore, synthetic data is constructed to demonstrate the 
preceding algorithm. 


As shown in Figure 2, the parent reflector has four separate areas of 
illumination on its surface. ^ Each area is assumed to have the same 
parabolic shape. For 0 0 £ ^ an< * e £ r < R 


r 


-^-r- [ k (H^) "2 (e) + ‘]- 


36 


u°<r,S) 


V* - r) 



(0) + i l» 


*]' 


17* 


q 3 ( e) + 1 , JT < 0 < 2 


130) 
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where 


( 0 ) * slnO + cos 0 


(31) 


The functions ^ 2 ( ® ) an< * qj(0) are cubic polynomial fits used to 

7T O <jy 

ensure smoothness in regions of 0 near 0 - — , *, — , 2*. Formulae for 

q 2 < 0 ) and 03 ( 0 ) may be found in [15] • The parameter k > 0, a 

stretch factor used to perturb the surface below the conic (k * 0 ) shape is 
taken as 0.25. 


For the complete surface, we define, for e £ r R, 


u(r,0) - ^ 


u^(r, 0 ) , 


■° ( r - e 'i) • 


0 , 


0 < 9 < - 


2 < © < ^ 


3* 


u (r ,0 - ir) , n < © < 


0 / 3ir\ 3w . „ ^ - 

^ <i ^r,0 - j-J, j~<Q< 2* 


(32) 


It is expected that the mesh will be stiffest near the outer hoop 
(r » R) and around the inner radius (r =* e), For this reason we choose a 
known value of E}(r) as 

E^r) - 2t - x sin |» ~ (e < r R) (33) 

* 

where t is a constant dependent on the mesh material. The stiffness in 
the angular d.lrecion is expected to oe uniform with 


E 2 (0) = x . (34) 

From data ptovid : 8 ] for the 100-meter point design, a reasonable 

value for t (given in units /N/m ) is 

t - 3.391 ; (35) 

similarly, other parameters are calculated to be uq ■ -7.5m, e » P,235m 
and R - 50m. 
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tu,(r£,0j) is calculated by 


A 10 x 24 grid of data points 
evaluating u(r,0) at points (r£,0j) with 

r i “ e + i — 7~ (1-1,2 L r - 10) (36) 


0 j = [7.5° + (j - >) 15° ] (j - 1,2 L e « 24) . (37) 

Values of 0j correspond to data taken along every other radial cord 
truss system with reflectors assumed located on the gore edges. 
Distributed loads are obtained by substituting (32)-(j4) int^ (1) and 
evaluating f(r ,0^. 

In the examples of the identification process to oe presented, an 
equal number of linear spline basis functions are used in both r end 0 
directions. That is, M * N + 1 for an increasing sequence of N-values. 
The cubic spline approximations (24) and (25) are ueed with fixed 
Ml “ Nj ■ 4 to represent Ej(r) and E?^). The IMSL version 
(ZXSSQ) of the Levenberg-Marquardt algorithm [17 j is employed to minimize 
jM>N given by (29). For the first choice of N ? nominal (v,6) parameter 
values to initialize the Levenberg-Marquardt algorithm are obtained by 
finding those (v ? 6) coordinates which cause (24) and (25) to best approxi- 
mate assumed functions E^(r) and E^(^) chosen as guessed forms for 
¥j(r) and Is^©), respectively. For larger N, the latest previously 

obtained set of converged coordinates is used as nominal parameters. 
Nr rical calculations are performed on a CDC Cyber 170-series digital 
;L.putv>r using default values of the IMSL convergence parameters. 

Two measures of Identification scheme performance are employed. The 
quantity 


} M ’ N - 



(38) 


is used as a measure o2 state estimati>n accuracy. Additionally, 


rM.N 



x 100Z 


(39) 
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measures the relative error between the true 


E( r ,0) - E 1 (r)E 2 (G) 
and the estimatec E(r,G) denoted by 

E M,N (r,0) - E“ ,N (r)E 2 * N (e) 


(40) 


(41) 


which is calculated from (24) and (25) using the (M,N) tli level of state 

approxit. rion obtained at step 8 of the computational procedure. In (39), 

i i M N 

y totes the norm on [e,R] x [0,2*]. R ’ provides a 

measure of parameter estimation accuracy. 


Convergence in the sense that 


rM.N 


-► 


0 


and 


?’ N 


-»• 


0 


as 


(M.N) * - 


depends on the ability of the cubic spline approximates (24) and (25) to 
accurately represent E^(r) and £ 2 ( 6 ). An exact pointwise fit can be 
obtained for ^(Q) by choice of the 4 6 -coefficients in (25). However, 
Ej(r) can at best be approximated to 


I E , ( r ) - E (r) | 

1.23Z 

E (r) 


relative error by (24) and (27) with M - 4, Consequently, entries in the 

M N N A 

(R ’ , J ’ ) sequence can be expected to cease decreasing past some (M.N) 

value. Less realistic examples in which i24) and (25) exactly fit simpler 


V 
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Ej(r) and £2(0) functions, and and R**’** monotonically decrease 

with increasing (M,N) can be found in [IS]. Also, using the best cubic 
spline fits to E^(r) and £^(9) obtained from (24) and (25) to 
define E(r,6), along with the exact f(r,G) data, we observed that 


?’ N - 0.087 


uniformly in (M,N). The following numerical results show that the 
parameter estimates froa the identification procedure tend to iaprove 

(reduce) this J M,N value at the expense of 

Example 1 : Estimate £3(0) boiling Ej(r) fixed at the best 

cubic spline estimate of E^fr) using (24). Nominal parameters for the 
N - 4 starting value are obtained by fitting (25) to 

E^O) - 1 + | cos 0 . 

Four 6-parameters are estimated and results summarized below. 


N 

?’ N , . 

r m ‘ n , z 

CP 

time, sec 

4 

0.0390 

5.13 


8 

6 

0.0384 

5.57 


23 

8 

0.0322 

5.69 


86 

10 

0.0347 

6.01 


105 

12 

0.0330 

5.83 


132 

Essentially no 

Improvement in state 

estimate was obtained 

past N - 8. The 

e M,N(q) ter< j e< j to 3*591 instead of 

E 2 (6 ) h 3.391 

• The • 0.20 bias Is 

attributed to 

the inability A (24) 

to exactly fit 

E^(r) 

• 
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Example 2 : Estimate Ej(r) holding E 2 (©) fixed at the best 
cubic spline estimate of using (25). Nominal parameters for the 
N * 4 starting value are obtained by fitting (24) to 


Four (►-parameters are estimated and results summarized below. 


N 

J M ’ N , m 

r m ‘ n , z 

CP time, sec 

4 

0.0355 

32.25 

22 

6 

0.0343 

24.5 

41 

8 

0.0270 

4.39 

75 

10 

0.0293 

13.17 

103 

12 

0.0275 

8.08 

130 

14 

0.0273 

7.44 

168 

16 

0.0271 

7.59 

222 

18 

0.0267 

7.68 

292 

20 

0.0264 

8.03 

370 

22 

0.0260 

7.91 

460 

24 

0.0267 

8.11 

578 

26 

0.0250 

7.49 

751 

28 

0.0203 

7.58 

847 

30 

0.0259 

7.71 

1050 

From 

a state estimation 

viewpoint, N ■ 

' 28 provides the best 

accuracy. 

Overall, considering 

state, parameter 

and ease of computation. 

N =■ 8 is 

values of 

best* Figure 3 shows the character of 
N* 

E^’^(r) for selected 
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Example 3 ; Estimate both Ei<r) and E 2 (0). Nominal parameters 

are obtained as before for N ■ 4 from 

Ej(r) = 5 

E^Q) - i - j sin 0 

For each N, the first coefficient, 6j, is held fixed at its initial 
value. Seven parameters are estimated. 


N 

J M ’ N , - 

R M * N , 2 

CP time, sec 

4 

0.0356 

32.24 

40 

6 

0.0341 

28.71 

67 

8 

0.0270 

4.42 

168 

10 

0.0293 

13.18 

209 

12 

0.0275 

8.09 

256 

14 

0.0273 

7.45 

337 

16 

0.0271 

7.59 

411 

IS 

0.0267 

7.69 

490 

20 

0.0264 

8.04 

567 

22 

0.0262 

7.90 

651 

24 

0.0260 

8.12 

768 

26 

0.0260 

7.47 

945 

Again, from 

overall considerations. 

N - 8 gives 

the best results 
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VII. CONCLUDING REMARKS 


In all examples we have been able to successfully estimate the surface 
shape of the model antenna. Simili" results have been obtained where 
random noise (approximately 5% noise level) has been added to the data. 
These and other findings may be found in Section VI of [15]. 
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Figure 1.- Side View of Maypole (Hoop/Column) Antenna. 



Figure 2.- Maypole (Hoop/Column) Antenna Reflector Surface. 
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Figure 3.- Estimate E^r) with E 2 (6) fixed. 
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